Giant optical nonlinearity of graphene in a strong magnetic field 
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We present quantum-mechanical density-matrix formalism for calculating the nonlinear optical 
response of magnetized graphene, valid for arbitrarily strong magnetic and optical fields. We show 
that magnetized graphene possesses by far the highest third-order optical nonlinearity among all 
known materials. The giant nonlinearity originates from unique electronic properties and selection 
rules near the Dirac point. As a result, even one monolayer of graphene gives rise to appreciable 
nonlinear frequency conversion efficiency for incident infrared radiation. 
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Graphene, a two-dimensional monolayer of carbon 
atoms arranged in a hexagonal lattice, holds many 
records as far as its mechanical, thermal, electrical, and 
optical properties are concerned; see. e.g. [T] for the 
review. With this Letter we would like to add yet an- 
other distinction to this list of superlatives: we show that 
graphene in a strong magnetic field has the highest in- 
frared optical nonlinearity of all materials we know. 

Strong optical nonlinearity of graphene, like most of its 
unique electrical and optical properties, stems from linear 
dispersion of carriers near the K,K' points of the Brillouin 
zone. As a result, the electron velocity induced by an in- 
cident electromagnetic wave is a nonlinear function of 
induced electron momentum. Nonlinear electromagnetic 
response of classical charges with linear energy disper- 
sion has been studied theoretically in [2] . Recently, four- 
wave mixing in mechanically exfoliated graphene flakes 
without magnetic field has been observed at near-infrared 
wavelengths [3]- Effective bulk third-order susceptibility 
was estimated to have a very large value, x ~ 1CP 7 
esu. This is comparable in magnitude to the resonant in- 
tersubband x^ nonlinearity observed in the mid-infrared 
range for low-doped quantum cascade laser structures [4] , 
which are essentially asymmetric coupled quantum well 
heterostructures. 

Nonlinear cyclotron resonance in graphene was consid- 
ered theoretically in [5], again in the classical limit, by 
solving the equation of motion F = dp/dt for a mass- 
less charge. Classical approximation can be applied to 
electrons in low magnetic field that are occupying highly 
excited Landau levels n> 1, when energy and momen- 
tum quantization are neglected. Here we present rigorous 
quantum mechanical description of the nonlinear opti- 
cal response of magnetized graphene, which is valid for 
arbitrary magnetic fields and electron distributions over 
Landau levels (LLs) . After finding matrix elements of the 
optical transitions between LLs, we calculate the third- 
order nonlinear susceptibility using the density-matrix 
formalism and then evaluate the efficiency of the four- 
wave mixing process. The magnitude of x turns out 



to be astonishingly large, of the order of 10 _1 esu at 
mid/far-infrared wavelengths in the field of several Tesla. 
This leads to a surprisingly high four-wave mixing effi- 
ciency of the order of 1CP 4 W/W 3 per monolayer. 

Linear (one-photon) absorption in monolayer and bi- 
layer graphene in arbitrary magnetic fields has been cal- 
culated in [6] using Keldysh Green's function formalism. 
This approach is inconvenient when it comes to calculat- 
ing the nonlinear optical response. The density matrix 
formalism adopted in this paper provides a rigorous, intu- 
itive, and straightforward framework for calculating the 
hierarchy of nonlinear optical susceptibilities and inter- 
action of strong multi-frequency EM fields or ultrashort 
pulses with graphene. Expressions for one-photon ab- 
sorbance obtained in [6] can be retrieved by calculating 
the linear susceptibility in the limit of a weak monochro- 
matic field. 

In the absence of the optical field, the effective-mass 
Hamiltonian [7HS] for a graphene monolayer (in the xy 
plane) in the magnetic field Bz, in the vicinity of K 
and K' points [TU] in the nearest-neighbor tight-binding 
model is given by 

/ TT X — ITT V \ 



\ TT X — ITTy / 

(1) 

where vp is a band parameter (10 s cm/s) [TTJ [T2], 7r = 
p + eA/c, p is the electron momentum operator, and A 
is the vector potential, which is equal to (0, Bx) here. 

To simplify notations, we write down the solutions to 
the Schrodinger equation Hq^ = e 1 ^ separately near the 
K and K' point. For example, near the K point the 
Hamiltonian is Hq = vf<j ■ if, where a = (a x ,a y ) is 
a vector of Pauli matrices. The cigenfunction is spec- 
ified by two quantum numbers, n and k y , where n = 
0, ±1, ±2, • • • , and k y is the electron wave vector along y 
direction [8]: 
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with C n = 1 when n — 0, C n = \j\[2 when n^O, and 



H H ({X - l 2 c ky)/l C ) 

il = / , , = — exp 



X Z' x I c ky 



where l c = \J ch/eB and H n (x) is the Hermite Polyno- 
mial. The eigen energy is e n = sgn(n)/kj c -y/jn|, where 
u c = V2v F /l c . 

In the presence of the incident classical optical field 
E = (l/2)eE ul e~ lut polarized along the vector e in the 
x-y plane (e LH s = [x - iy]/y/2 and e RHS = [x + iy]/V%, 
which denote left-hand and right-hand circularly polar- 
ized light) , we add the vector potential of incident optical 
field, A opt = ~E, to the vector potential of the mag- 
netic field in the generalized momentum operator if in 
the Hamiltonian. This results in adding the interaction 
Hamiltonian Hi nt to H$, where 



c 



(3) 



This linear in A opt expression for the interaction 
Hamiltonian is exact, unlike the case of the kinetic en- 
ergy operator quadratic in momentum, where the term 
proportional to A 2 is usually neglected. Note also that 
Eq. ([3]) does not contain the momentum operator and 
its matrix elements are simply determined by the matrix 
elements of a. This immediately gives the selection rules 
[6] for the transitions between the LLs: clhs photons are 
absorbed when \nf \ = |nj| + l, whereas Hrhs photons are 



absorbed when 



l f\ 



1. Here n,- and 



indicate 



initial and final energy quantum numbers of LLs. 

Now we can write a standard time-evolution equation 
for the density matrix of Dirac electrons in graphene cou- 
pled to an arbitrary optical field: 



^ = -^[H + H inU p]+R(p). 



(4) 



Here R(p) describes incoherent relaxation due to disor- 
der, interaction with phonons, and many-body carrier- 
carrier interactions. Equations Eq. Q have to be solved 
together with Maxwell's equations that contain the opti- 
cal polarization P(r,t) = (l/V)Tr(/5 • p) (average dipole 
moment (ft) per unit volume) as a source term. In the 
perturbative regime, they give rise to the hierarchy of the 
optical susceptibilities x [T3] . but they are also valid 
for describing non-perturbative coupling to strong fields, 
interaction with ultrashort pulses etc. 

Since graphene is essentially a 2D system, it makes 
sense to introduce a surface (2D) polarization P s deter- 
mined as an average dipole moment per unit area rather 
than unit volume. Below we will use 2D susceptibilities 
unless specified otherwise. 

For a weak monochromatic field one can retain only 
the term p\ 



(i) 
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J2m n Pnmfimn to obtain an expression for the linear sus- 
ceptibility: 
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ifolmn) linear with respect to the field and take the sum 



Here we used (m\Hi nt \n) = —(i/u})evF(m\a\n)E(uj) 
and (m\jl\n) = (ih/(e n — e m ))evp (m\a\n). Note that 
the matrix element of the interaction Hamiltonian can 
be written as —p mn E, where p mn = (i/co)evF(iTi\(f\n) , 
and fx mn = p mn when e n - e m — frw. 

We assumed for simplicity that the relaxation term 
for the off-diagonal density matrix elements R m n = 
—ImnPmn and all 7's are the same. For easy comparison, 
we used the same notations for LLs as in [5]: a, d = ± 
denote whether the corresponding state belongs to the 
conduction (+) or valence (-) band and v n ,a are the fill- 
ing factors of LLs; a complete occupation corresponds to 
v = 2. The degeneracy of a given LL is 2/(tt1 2 ) including 
both spin and valley degeneracy. After calculating the di- 
mensionless linear absorbance as (27rcj/c)Im[x^ 1 - ) (w)] we 
obtain the same result as in [B]. 

Now we consider a specific example of the nonlinear op- 
tical interaction, namely the four-wave mixing. Consider 
a strong bichromatic field E = (1/2)(E\ cxp(— iuit) + 
E2 exp(— ioj%t) + c.c.) normally incident on the graphene 
layer. Here u)\ is nearly resonant with the transition from 
n = — 1 to n = 2 and E\ has left circular polarization. 
The frequency u>2 is nearly resonant with the transition 
from n = to n = ±1 and E2 has linear polarization, 
so that it couples both to transition — 1 — ► and — s- 1, 
as shown in Fig. 1. As a result of the four- wave mixing 
interaction, the right-circularly polarized field E3 at fre- 
quency UJ3 — u>i — 2ui2 nearly resonant with the transition 
from n = 2 to n — 1 is generated. 

Efficient nonlinear mixing becomes possible due to 
strong non-equidistancy of the LLs and unique selection 
rules A|n| = ±1 which enable transitions with change in 
n greater than 1, for example the transition from state 
n = — 1 to state n — 2. This transition would be for- 
bidden in conventional LL systems with An = ±1 selec- 
tion rule. The effective dipole moments for all transitions 
shown in Fig. lb scale as np/w, i.e. they are similar to 
each other within a factor of 2 and are very large: of 
the order of 10-100 A in the mid/far-IR range. This, in 
combination with sharp peaks in the density of states at 
LLs enables a strong nonlinear response. 

For simplicity we assume that the incident field is 
not strong enough to significantly modify populations 
and all states below n = are fully occupied. Then 
the optical fields interact resonantly only with states 
n = —1, 0, 1, 2, which we renamed to n = 1, 2, 3, 4 in Fig. 
lb. The Hamiltonian can be truncated to a 4x4 matrix, 
where (Ho) mn is diagonal, with diagonal elements being 
the energies of corresponding LLs, and the interaction 
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FIG. 1: (a): Landau levels near the K point superimposed on 
the electron dispersion without magnetic field E — ±vp|p|. 
(b): A scheme of resonant four- wave mixing process in the 
four-level system of LLs with energy quantum numbers n = 
— 1,0, +1,4-2. The case of exact resonance is shown. Polar- 
ization of light corresponds to the allowed transitions. 



Hamiltonian is given by the matrix —\x mn E as specified 
above. This approximation is similar to the one adopted 
in ^ I14H16) for analyzing resonant nonlinear processes 
in coupled quantum-well heterostructures. The resulting 
third-order nonlinear optical susceptibility at frequency 
CJ3 = uj\ — 2w2 is 



x (3) M 
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P33 

r 



(^) 3 r 43 
P22 , P22 — pn 



31 1 32 
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(6) 



Here the complex detuning factors are T2i = 721 +i{{s2 — 
ei)/h - w 2 ), T 32 = 732 + «((£3 - - w 2 ), r 4i = 
74i+i((e4-£i)/^-wi), and T 43 = 743-H((e 4 -£3) /h-(J 3 )- 
In deriving Eq. Q from Eq. (|4 ) we assumed that pop- 
ulations Pram are constant and solved for the off-diagonal 
density matrix elements. To get an order-of-magnitude 
estimate, we assume that all fields are in exact resonance 
and all dephasing rates are the same, so that = 7. 
We also assume for definiteness that state 1 is fully oc- 
cupied while states 2,3, and 4 are empty. Coming back 
to original notations of LLs in Fig. la, this means that 
the n — LL is empty, i.e. the Fermi level is between 
states n = and n = — 1. The magnitude of x' 3 - 1 will 
be similar for any distribution of carriers as long as not 
all population differences in Eq. (|6| are zero. Then we 
obtain 



(3) 
X26 



(6/7r^)/i 43 /i4i//32/221 



3.6x10" 



1 



B(T) 



(7) 



where the magnetic field B(T) is expressed in Tesla and 
we took 7 = 2 x 10 13 s _1 for the dephasing rates. 



This is a 2D (surface) susceptibility. To compare 
with bulk materials, we divide Eq. |7| by the mono- 
layer thickness ~ 3 A to obtain the bulk susceptibility 
X3D ~ 1 {^/B(T)) esu. This is by far the strongest non- 
linearity as compared to any material that we are aware 
of. The frequencies involved in the four- wave mixing pro- 
cess fall into the mid/far-IR range for the magnetic field 
of a few Tesla, as shown in Fig. 2. In particular, at B = 1 
T the generated nonlinear signal is at the wavelength of 
about 82 fim. 

The magnitude of the x^ 3 "* nonlincarity scales roughly 
as I/7 3 , i.e. it rapidly decreases with increasing line 
broadenings. However, even in a very disordered mate- 
rial with broadenings of the order of transition frequen- 
cies ~ 10 14 s _1 the magnitude of is still record-high: 
above 10~ 5 esu. 
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FIG. 2: Transition frequencies in the above 4-energy level 
graphene system. Wjj indicates the transition frequency be- 
tween level i and j. 

From wave equations, the electric field amplitude 
E^uis) of the generated nonlinear signal is given by 

E 3 = i(27rw3/c)x 2 3 D(w3)Si£ ; 2- Assuming that all beams 
ideally overlap within the area A of the graphene sample, 
the power of the nonlinear signal from one monolayer is 



P(w 3 ) 



167r 2 a; 3 
c 2 A 



(x (3) ) 2 P 1 K)(P 2 (^ 2 )) 2 . (8) 



For the illuminated area A = 10 4 cm 2 , the power 
conversion efficiency for the nonlinear signal generation 
scales as - 10" 4 (1/B(T)) W/W 3 . This is a remarkably 
large efficiency for one monolayer of material. It can be 
further increased by stacking several layers of graphene, 
e.g. by fabricating non-Bernal stacked epitaxial graphene 
layers as in recent demonstration of light amplification in 
graphene [T7] , 

The above expression for the nonlinear power becomes 
invalid at very high optical fields, when the effective Rabi 
frequencies fl mn E/H of the optical fields become of the or- 
der of y/^/Tmn, where T mn is population relaxation time 
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between states m and n. For the magnetic field of sev- 
eral Tesla and r ~ 1 ps, this corresponds to intensities 
of about 10 4 W/cm 2 . At higher incident fields, popu- 
lation differences in Eq. ^ become optically saturated 
and eventually decrease as 1/E 2 . As a result, the growth 
of the nonlinear power with increasing incident power 
slows down. At even higher optical fields, when the Rabi 
frequencies exceed the line broadenings 7, the broaden- 
ing factors T's at the transitions driven by strong optical 
fields increase as E 2 and the nonlinear power decreases 
with increasing incident power Pi ~ P 2 ~ P as 1/-P- To 
analyze the strong-field case quantitatively, one needs to 
solve Eq. Q for both diagonal and off-diagonal density 
matrix elements, which is more tedious but straightfor- 



ward. 

In conclusion, graphene in a strong magnetic field 
possesses record-high optical nonlinearity due to unique 
properties of quantized Landau levels near the Dirac 
point and selection rules for the optical transitions be- 
tween Landau levels. High nonlinearity leads to signif- 
icant nonlinear frequency conversion efficiency even for 
one monolayer of material. The nonlinearity is expected 
to be ultrafast, enabling response to THz modulation. 
These unique properties of magnetized graphene may 
have important implications for coherent nonlinear gen- 
eration and detection in the mid-infrared and THz range. 
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